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Abstract 

We study x c i(2P) inclusive hadroproduction at next-to-leading order (NLO) within the fac- 
torization formalism of nonrelativistic quantum chromodynamics (NRQCD), including both the 
color-singlet 3 p| and color-octet 3 jS{ cc Fock states. Assuming the recently discovered X(3872) 
hadron to be the 2P (1 ++ ) charmonium state, we perform a fit to the cross sections measured by 
the CDF, CMS, and LHCb Collaborations. We either obtain an unacceptably high value of \ 2 or 
a value of I-R^pWI incompatible with well-established potential models. We thus conclude that 
NLO NRQCD is incompatible with the hypothesis X(3872) = Xci(2P). 

PACS numbers: 12.38.Bx, 12.39. St, 13.85.Ni, 14.40.Pq 
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During the past decade, a series of charmonium or charmonium-like X, Y, Z states were 
discovered (for a recent review, see Ref. jj). The X(3872) state is one of the most interesting 
among them. It was discovered in 2003 by the Belle Collaboration at KEKB in B meson 
decays {2], and confirmed shortly afterwards by the BaBar Collaboration at SLAC PEP-II js|. 
It was also observed by the CDF [4] and DO [5] Collaborations in pp collisions at the Tevatron 
Fermilab. Ever since its discovery, many theoretical group have tried to interpret its nature, 
which has remained mysterious to date, and it is an urgent task of great importance and 
broad interest to solve this notorious puzzle of hadron spectroscopy. Typical options include 
conventional charmonia [6], D*°D° /D°D*° molecules and tetraquarks [S|. However, none 
of them can provide a convincing description of all the experimental measurements. After 
analyzing the dipion mass spectrum in X73872) — > J /if) + vr + 7r~, only two options for its 
jpc p r0 p er ty are left, either 1 ++ or 2~ + [9]. a In pp and pp collisions, most of the X(3872) 



mesons are produced promptly rather than through decays of b hadrons [ll|,[12|]. The study of 
X(3872) prompt production provides complementary information on its nature. In Ref. 13], 
the cross section of X(3872) was estimated under the assumption that it is a loosely-bound 
D*°D /D°D molecule, and the upper bound of the theoretical calculation was found to 



be much smaller than the CDF measurement 



111 ] . Later, Artoisenet and Braaten 



pointed out that the upper bound of this prediction can be rendered consistent with the 



flfl 



ll| by properly taking into account rescattering effects. They also used 
actorization approach [15| to interpret X(3872) prompt 
11| and presented predictions for the LHC. However, their 
oredictions significantly exceed the new measurements reported by the CMS 1^] and LHCb 



Tevatron data 
the nonrelativistic QCD (NRQCD) 

nn 

production at the Tevatron [4J, [5|, 



16[ Collaborations. In their charmonium interpretation [14J, X(3872) is assumed to be a 
1 ++ state that is dominantly produced via the color-octet cc Fock state 3 Sf \ and the short- 
distance coefficients are calculated at leading order (LO). Note that, at first sight, the mass 
value 3.872 GeV seems too low for a x cl (2P) candidate, but color- screening effects together 
with coupled- channel effects may draw its mass down towards 3.872 GeV [171 ]. 

Recent NRQCD analyses have revealed that next-to-leading order (NLO) corrections play 
a key role in explaining the J /if} 18] and % c j(lP) [19( yields measured at the Tevatron and 



Very recently, the LHCb Collaboration [l0( established the assignment J PC — 1 ++ , which, however, still 
lacks independent confirmation. 
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FIG. 1: Overview of the IR singularity structure. 

the LHC. Under the assumption that X(3872) is the 1 ++ charmonium state x cl (2P), it is 
then natural to ask if its prompt production rates may be explained upon including NLO 
corrections. The main goal of our work is to answer this question. To this end, we shall first 
calculate the cross section of inclusive x c i(2P) hadroproduction at NLO in NRQCD and 
then check if its free parameters can be adjusted so as to yield a satisfactory description of 
the available prompt X(3872) hadroproduction data 



Invoking the factorization theorems of the QCD parton model and NRQCD [l5] 
inclusive x cl (2P) hadroproduction cross section is evaluated from 

da(AB -> X ci(2P) +X) =YlJ dxdyf t/A (x)f j/B (y) 



the 



x (0 Xcl{2P) [n})da{ij ->■ cc[n] + X), (1) 



where fi/A( x ) are the parton distribution functions (PDFs) of hadron A, {0 Xcl ^ 2P \n\) are 
the long-distance matrix elements (LDMEs), and da(ij — > cc[n] + X) are the partonic cross 
sections. Working in the fixed-flavor-number scheme, % and j run over the gluon g and 
the light quarks q = u,d,s and anti-quarks q. The system X always contains one hard 
parton at LO and is taken to be devoid of heavy flavors, which may be tagged and vetoed 
experimentally. The contribution due to final states in which X comprises an open cc pair 



is found to be suppressed by one order of magnitude [20]. At LO in the relative velocity v 
of the bound c and c quarks in the charmonium rest frame, only the states n = 3 p\ , 3 Sf 
contribute, the contributions due to other states being at least of relative order 0(v 4 ) |15| . 

In our calculation all singularities are canceled analytically. The ultraviolet divergences 
are removed by renormalizing the parameters a s and m c and the wave functions of the exter- 
nal lines. Specifically, we work in the on-shell scheme, except for a s , which is treated in the 
modified minimal-subtraction scheme. The infrared (IR) singularities are canceled similarly 



as described in Ref. 



21] . Fig. [TJ displays an overview of the IR singularity structure. In the 



case of the 3 Sj 8 ' state, the soft and collinear singularities of the real corrections are canceled 
by complementary contributions stemming from the virtual corrections up to universal parts 
to be absorbed into the hadron PDFs. In case of 3 P[^ state, the real corrections produce 
two new classes of soft singularities, named soft #2 and soft #3 terms, on top of the soft #1 
terms familiar from the 3 S± channel. The soft #2 terms are not proportional to LO cross 
sections; they cancel against the IR singularities of the virtual corrections left over upon 
the familiar cancellation against the soft #1 terms. The soft #3 terms cancel against IR 
singularities stemming from radiative corrections to the (0 Xcl ^ 2P \ 3 Si )) and {0 Xcl ^ 2P \ 3 Sf^)) 
LDMEs. Notice that the inclusion of the 3 Sf^ contribution is indispensable in order to obtain 
an IR finite result. We thus recover the notion that the color-singlet model is not a complete 
theory. By the same token, the dependencies of (0 Xcl ^ 2P \n\) and da(ij — > cc[n] + X) on the 
NRQCD factorization scale fi\ only cancel after summation over n. 

We now describe the choices of input for our NLO NRQCD calculation. We take the 
charm-quark mass to be m c = 1.5 GeV and use the two- loop formula for a^ 1 ^ with n/ = 4 



22 1, which comes 



active quark flavors. As for the proton PDFs, we adopt the CTEQ6M set 
with asymptotic scale parameter Aq^ d = 326 MeV. We choose the MS renormalization, 
factorization, and NRQCD scales to be fi r = /if = ^tut and /ia = V m c, where itlt = 
\J Px + 4m;? is the Xci(2-P) transverse mass, and independently vary £ and r/ by a factor of 
two up and down about their default values £ = rj = 1 to estimate the scale uncertainty. To 



LO in v, we have 



15| 



(OXe 1 (2P) ( 3p[il )) = ( 2j + 1 )^\ H ^ p (o)^ ( 2 ) 

where Ca = N c = 3, J = 1 is the total angular momentum of the % c i(2P) meson and B,2p(r) 
is its radial wave function, which may be calculated using models for the QCD potential 
of the charm quark. Adopting frequently used potential models with different choices of 
parameters, |i^ P (0) | 2 is found to range from 0.076 GeV 5 to 0.183 GeV 5 Q. As the default 
for our fits, we adopt the value l-R^f (0) I 2 = 0.102 GeV 5 obtained using the Buchmiiller- 
Tye potential 24j]. To compare theoretical predictions with the experimental data, we also 
need to know the branching fraction (BR) of the decay mode X(3872) — > J/ip + 7r + 7r~ 
used to identify the V(3872) meson. It has not been determined yet, but the lower bound 



BR > 2.6% has been established at 90% C.L. |25j. Furthermore, the upper bound BR < 9.3% 

n 

was derived at 90% C.L. using constrains from some other decay channels [14] . In our fits, 



we use BR = 2.6%. 

Based on the measurements by the CDF Collaboration 



111 ], at center-of-mass energy 



\/s = 1-96 TeV, the prompt production cross section of X(3872) mesons with rapidity 
\y\ < 0.6 and transverse momentum px > 5 GeV is estimated to be \v\ 

a%™ p \pp -> X(3872) + X)BR = (3.1 ± 0.7) nb. (3) 

At LHC, prompt X(3872) production was first measured by the CMS Collaboration [12j], at 
y/s = 7 TeV, with the result 

ag r Ms Pt (pp -> X(3872) + X)BR = (1.06 ± 0.19) nb, (4) 



for \y\ < 1.2 and 10 GeV < pt < 30 GeV. They also presented a pt distribution [12 ]. 



The LHCb Collaboration also measured X(3872) production at yfs = 7 TeV, but did not 
discriminate between 6-hadron and prompt sources [16]. Averaging the CDF and CMS 
measurements of the non-prompt fraction, (16.1 ±4.9±2.0)% j^, 11] and (26.3±2.3±1. 







12j, respectively, we estimate the LHCb prompt cross section to be 

<H°cb Pt (M° -> ^(3872) + X)BR = (4.26 ± 1.23) nb, (5) 

for 2.5 < \y\ < 4.5 and 5 GeV < p T < 20 GeV. 

In the following, we perform a NLO NRQCD test of the hypothesis that the X(3872) 
hadron is the x cl (2P) charmonium state. Since V(3872) production via feed-down of heavier 
charmonia has not been observed, we assume for the time being prompt production to be 
approximately exhausted by direct production. In fact, charmonia heavier than the V(3872) 
hadron have sufficient phase space above the open charm production threshold and preferably 
decay to pairs of D mesons so as to evade the kinematic constraint of cc bound state 
formation. We are thus led to fit Eq. (CQ), which depends on the parameters (0 Xcl ^ 2P H 3 p| 1 ^)) 
and (e> Xcl(2P) ( 3 Sf ] )), to the experimental data of prompt V(3872) production |j, [ll|, LL2, M • 
We consider four options altogether. On the theoretical side, we either fix (0 Xcl ^ 2P \ 3 P^)} = 
0.438 GeV 5 by the potential model of Ref. |24j, or fit it along with (0^ l( - 2P \ 3 Sf ] )) . In the 
latter case, we actually take (O x ^ 2P \ 3 P[ 1] ))BR and (O x ^ 2P \ 3 Sf ] ))BR to be fit parameters. 
On the experimental side, we either include the LHCb result [16[ of Eq. (jSJ) in the fit along 
with the CDF result 0, LLlJ of Eq. ((3j) and the CMS measurement of the pt distribution 
12J, which includes four data points, or we exclude it. In order to avoid double counting, 
we always exclude the CMS result [12j of Eq. (@| from the fit. 



TABLE I: Results of our NLO NRQCD fits to the measured p T distribution of prompt X(3872) 
production [2] and the integrated cross section of Eq. © including or excluding the result of 
Eq. ©. In the one-parameter case, (0 Xcl ^ 2P \ 3 P^)) is determined by the potential model of 
Ref. [24]], while it is a fit parameter in the two-parameter case. We adopt BR = 2.6%. 





One-parameter case 


Two-parameter case 




w/ LHCb data 


w/o LHCb data 


w/ LHCb data 


w/o LHCb data 


(0Xcl(2P)(3p[l])) [Q eV 5] 


0.438 


0.438 


n in+0.050 
u - iu -0.050 


n iq+0.092 
u - ly -0.094 


(e^( 2p )( 3 sf)) [GeV 3 ] 


(3.84±8;||) x 10~ 3 


(4.30l° : 3 °) x 10- 3 


(2.95l°;i) x 10- 3 


(3.36ig:|g) x io- 3 


a—BR [nb] 


4.86^ : g 


5.06ll;l 


2 05 +a33 

Z - UO -0.33 


o qo+0.48 


^MS Pt BR [nb] 


1.26l°; 2 ° 




1.041^ 


i-i6±g:S 




1 3 fifi +317 
10.DD_2 .65 


14.2ll|g 


5.69l°jJ 


8.04ii:l 


X 2 /d.o.f. 


79.1/5 = 15.8 


16.7/4 = 4.18 


4.26/4 = 1.07 


0.63/3 = 0.21 



The results of the four fits are summarized in Table [H and their goodness may be conve- 
niently assessed from Fig. |2j In the two-parameter case, the quoted values of (£> Xcl(2P) ^Pf 1 )) 
and (0 Xcl ( 2P \ 3 Sf^)) correspond to our default value BR = 2.6%. When Eq. ([5]) is excluded 
from the fit, dLy^^BR is a genuine prediction. The uncertainties are estimated by adding 
in quadrature the errors of experimental origin resulting from the fits using our default NLO 
NRQCD results and those due the variations of the scale parameters £ and r\. 

The one-parameter fit including the LHCb data point of Eq. ((5]) has d.o.f . = 5 degrees of 
freedom and yields x 2 — 79.1, so that xVd.o.f. = 15.8 is intolerably large suggesting that the 
experimental data is poorly described by only adjusting {0 Xcl ^ 2P ^ ( 3 S^)) . This is also evident 
from the upper left panel in Fig. [2j We observe that x 2 rapidly increases with |-R'2p(0)| 2 
and BR. For |i2^ P (0)| = 0.076 GeV 5 and BR = 2.6%, we obtain the best value x 2 = 37.5, 
which is still unacceptably large, while for |-R 2 p(0)| 2 = 0.183 GeV 5 and BR = 9.3%, x 2 is 
around 4900. We also notice that the NLO NRQCD result greatly overshoots the LHCb 
data point although it is included in the fit. Excluding it from the fit mildly increases the 
central values of the fit results and their errors, but significantly reduces x 2 1 by almost a 
factor of five. The NLO NRQCD prediction for Eq. © is then nb, so that the 

theory band overshoots the LHCb result by almost six experimental standard deviations. 

The two-parameter fit including the LHCb data point of Eq. (J5J) works nicely, yielding 
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xVd-o.f. = 1.07. Specifically, we have (0 Xcl ( 2P) ( 3 P[^ ) ) BR = (2.601^) x 1(T 3 GeV 5 and 



(0^ 2P \ 3 Sf ] ))BR = (7.67±}f ) x 10- 5 GeV 3 . From BR > 2.6% at 90% C.L. 25], we thus 
derive the 90% C.L. upper bound \R' 2P {0)\ 2 < (2.33± 1.16) x 10~ 2 GeV 5 , which undershoots 
the smallest known potential model result, 0.076 GeV 5 [23], by more than 4a and whose 
central value is more than a factor of three smaller than the latter. Conversely, if we choose 
|i?2p(0)| 2 within the ballpark of potential model calculations, then the upper bound on BR 

I L 

is around three times smaller than the lower bound 2.6% [25]. Also the LHCb data point is 
nicely described by the fit. Excluding it from the fit appreciably increases the central values 
and errors of the fit results and pushes xVd.o.f. far below unity, to 0.21. Specifically, we 
find (C^( 2P )( 3 PJ 1] ))BR = (4.94±^) x lO" 3 GeV 5 and (0^ 2P \ 3 S f ] ))BR = (8.74+^) x 
10~ 5 GeV 3 . However, the NLO NRQCD prediction for Eq. © now reads (8.04+};^) nb, 
so that the theory band overshoots the LHCb result by almost two experimental standard 
deviations. 

In conclusion, we tested the hypothesis that the X(3872) hadron, whose nature is re- 
maining undetermined even a decade after its discovery, is a pure Xci(2-P) charmonium 
state, by^ fitting all available data of prompt AT(3872) production, from the CDF {4], 11]. 



CMS 12j, and LHCb [16| Collaborations, at NLO in the effective quantum field theory 



of NRQCD endowed with the factorization theorem proposed by Braaten, Bodwin, and 
Lepage [ijj. NRQCD factorization, which is arguably the only game in town among the 
candidate theories of heavy-quarkonium production and decay, has recently been impres- 
sively consolidated at NLO by global analyses of the world data of J/ip inclusive production 
(for a review, see Ref. [26i]). Assuming the color-singlet LDME (0 Xcl( - 2P \ 3 Pl^)) to be in 
the ballpark of well-established potential models 23[ and imposing the lower bound on 
the BR of X(3872) — > J/ip + 7r + 7r~ quoted by the Particle Data Group 25j, we find that 
the pure x cl (2P) assignment to the V(3872) hadron is strongly disfavored. The tension 
may be somewhat relaxed by excluding the LHCb data point 16] from the fit, which is, 



however, unmotivated and unsatisfactory, the more so as this challenges the CDF [J, [H] 
and CMS [3] measurements of the non-prompt V(3872) BR. In order to obtain a good 
fit, one is forced to sacrifice Eq. (J2]), one of the tenets of NRQCD factorization 15|] , and 
to allow for (0 Xcl ( 2P )( 3 p| 1 ')) to be an independent fit parameter. The price to pay is that 
comes out more than a factor of three, which amounts to more than 4<r, 



below any known potential model result |23j in combination with the BR lower bound 



25]. 
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FIG. 2: The prompt A(3872) production cross sections measured by the CDF 
and LHCb [16J Collaborations are compared with NLO NRQCD results based on one-parameter 
(upper row) or two-parameter (lower row) fits including (left column) or excluding (right column) 
the LHCb data point of Eq. ([5]) [16]. The dotted, dashed, and solid lines represent the 3 pj^ and 
contributions and their sum, respectively. The yellow bands indicate the uncertainties in the 
total results. 



Despite concerted experimental and theoretical endeavors during the past decade, the quest 
for the ultimate classification of the X(3872) resonance remains one of the most tantalizing 
puzzles of hadron spectroscopy at the present time. 

This work was supported in part by the German Federal Ministry for Education and 
Research BMBF through Grant No. 05H12GUE. 
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